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ABSTRACT 



We present the results of a series of XMM-Newton EPIC and OM observations 
^ • of Her X-l, spread over a wide range of the 35 day precession period. We confirm 

00 I that the spin modulation of the neutron star is weak or absent in the low state - in 

0^ ■ marked contrast to the main or short-on states. During the states of higher intensity, 

we observe a substructure in the broad soft X-ray modulation below ~ 1 keV, revealing 
the presence of separate peaks which reflect the structure seen at higher energies. The 
strong fluorescence emission line at ~6.4 keV is detected in all observations (apart from 
one taken in the middle of eclipse), with higher line energy, width and normalisation 
during the main-on state. In addition, we report the detection of a second line near 
^ ■ 7 keV in 10 of the 15 observations taken during the low-intensity states of the system. 

Oh| This feature is rather weak and not significantly detected during the main-on state, 

when the strong continuum emission dominates the X-ray spectrum. Spin resolved 
£^ ■ spectroscopy just after the rise to the main-on state shows that the variation of the Fe 

^ \ Ka at 6.4 keV is correlated with the soft X-ray emission. This confirms our past finding 

based on the XMM-Newton observations made further into the main-on state, and 
indicates the common origin for the thermal component and the Fe Ka line detected 
at these phases. We also find that the normalisation of the 6.4keV line during the low 
state is correlated with the binary orbital phase, having a broad maximum centered 
near 4> rbit ~ 0.5. We discuss these observations in the context of previous observations, 
investigate the origin of the soft and hard X-rays and consider the emission site of the 
6.4keV and 7keV emission lines. 

Key words: accretion, accretion disks - X-rays: binaries - individual: Her X-l - 
stars: neutron 



1 INTRODUCTION 



Her X-l is one of the best studied X-ray binaries in the 
sky; this is borne out by the fact that since 2000 there has 
been more than 30 papers published in refereed journals 
with "Her X-l" in the title. We list very briefly the most 
salient observational characteristics which are relevant to the 
present study. Her X-l is a binary system, which consists of 
a neutron star and an A/F secondary star; it has an orbital 
period of 1.7 day and the spin period of the neutron star 
is ~1.24 sec. It is one of only a few systems which shows a 
regular variation in X-rays: the 35 day variation is referred to 



as the "beat" period throughout this paper, and the phasing 
on this period is marked as $35. 

The beat cycle begins with a sharp rise in X-ray in- 
tensity, which reach the maximum flux, F max , in less than 
~0.5 day and then lasts ~10 day ("main-on" state). This is 
followed by a period of low X-ray emission ( "low-state" ) and 
later again by a "short-on" phase, during which the X-ray 
flux reaches ~ l/20F ma x and ~ l/SF max , respectively. A 
second low-state follows until the onset of the next main-on. 
The spin modulation of the neutron star is most prominent 
during the main-on state, becoming much less so (if seen 
at all) during the low-state and increasing again during the 
short-on (see e.g. Deeter et al. 1998). 
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The observed X-ray spectrum changes markedly over 
the beat cycle. However, it has been recognized that the X- 
ray spectra taken at different beat phases can be fitted using 
a constant emission model and adding one or more absorp- 
tion components which vary in column density over the beat 
period, reaching a maximum value during the low-state (e.g. 
Mihara et al. 1991, Leahy 2001). This strengthened the idea 
that the 35 day cycle is due to the precession of an accretion 
disk, that periodically obscures the neutron star beam. 

A rather broad Fe Ka emission line at ~6.5 keV has 
been observed by Ginga (see Leahy, 2001 and past references 
therein); line energy and width have been found to vary over 
both the beat and the spin period (e.g. Choi et al. 1994). 
Using ASCA data taken during the main-on state, Endo, Na- 
gase & Mihara (2000) were able to resolve the feature into 
two discrete emission lines, at ~6.4 keV and 6.7 keV. How- 
ever, the second feature has never been detected in other 
observations. With the advent of imaging X-ray satellites 
with large effective area over a broad energy band and rel- 
atively good spectral resolution at energies near 6 keV, it 
is now possible to resolve the Fe lines with better signal to 
noise than what was possible using ASCA. Moreover, while 
ASCA had a lower limit of ~0.5 keV, the energy response of 
XMM-Newton extends down to ~0.2 keV and this allows a 
more accurate quantification of the absorption column den- 
sity. 

Another manifestation of the 35 day cycle is the evolu- 
tion of the Her X-l pulse profile at energies above ~ 1 keV. 
A few other pulsars show pulse shape changes which are cor- 
related with high-low states of X-ray intensity, but they are 
typically transients such as GX 1+4 (Dotani et al. 1984) and 
EXO 2030+375 (Parmar, White & Stella 1989) or flaring 
sources as LMC X-4 (Levine et al. 1991). Her X-l is unique 
in having a repetitive evolutionary pattern in pulse changes 
tied to a regular high-low intensity cycle, and is therefore a 
subject of continued interest. Particularly successful in ex- 
plaining the basic features observed by Ginga and RXTE 
is the accretion column model by Scott et al. (2000). This 
model invokes the successive obscuration of a compact pen- 
cil beam from the neutron star, two fan beams with size 
~ 2R ns and a more extended scattered component. 

At lower energies (below ~ 0.7 keV), the light curve 
does not show evidence for a similar variety of emission 
components. Past detections of pulsations in the soft X- 
rays of Her X-l obtained using EXOSAT (Ogelmen & 
Triimper 1988), ROSAT (Mavromatakis 1993), and Bep- 
poSax (Oosterbroek et al. 1997, 2000, 2001) only show a 
broad and quasi-sinusoidal modulation, which is shifted in 
phase with respect to the maxima observed at higher ener- 
gies. The phase shift, as measured during the main-on and 
short-on state, differs from 180° . Under the assumption that 
soft X-rays are due to the reprocessing of the pulsar beam 
by the inner edge of the disk, this is usually interpreted as 
evidence for a tilt angle in the disk (see e.g. Oosterbroek 
at al. 1997, 2000). As predicted by precessing disk mod- 
els (Gerend & Boynton, 1976), the hard/soft shift in phase 
should vary along the beat cycle. The ideal way to test this 
would be to track the pulse difference over the entire 35 day 
period. However, past measurements were restricted to the 
main-on and short-on states, due to the difficulty of detect- 
ing a spin modulation during phases in which the system 
is highly absorbed. Using XMM-Newton data, Ramsay et 



al. (2002, hereafter R02) reported the first evidence for a 
substantial change in the phase difference along the beat 
cycle. We also found that the values of the phase shift ob- 
served by XMM-Newton during the main-on and short-on 
differ from those measured in the past at similar beat phases. 
However, at the time of that paper only three XMM-Newton 
observations, at three different beat phases, had been per- 
formed. 

To date, Her X-l has been observed by XMM-Newton 
on 15 separate occasions giving good coverage over the beat 
period. The first three observations have been discussed in 
detail by R02 and Jimenez-Garate et al. (2002). Here we 
report the analysis of the remaining datasets, focusing on 
results obtained using the broad energy band EPIC detec- 
tors and putting our past findings in a broader contest. We 
also present the simultaneous UV observations made using 
the Optical Monitor (OM). Observations are described in 
§ 2. UV data and timing analysis are presented in § 3 and 
§ 4 respectively, while the behaviour of the Fe feature(s) over 
the beat, the spin and the ornital period is reported in § 5. 
Summary and discussion follow in § 6. 



2 OBSERVATIONS 

2.1 Instrumental setup and data reduction 

XMM-Newton has three broad-band imaging detectors with 
moderate spectral resolution: two EPIC MOS detectors 
(Turner et al. 2001) and one EPIC PN detector (Striider 
et al. 2001). During the observations of Her X-l, PN and 
MOS1 were configurated in timing mode, while MOS2 was 
configured in either small window or full frame mode. All 
three EPIC cameras were used in conjunction with medium 
filters. The Optical Monitor (Mason et al. 2001) was con- 
figured in imaging mode, by using the UVW1 and UVW2 
filters (effective wavelength 2910 and 2120 A respectively). 

The details of the observations are summarized in Ta- 
ble 1; as previously mentioned, the results from the first 
three observations have been already reported in detail by 
R02 and Jimenez-Garate et al. (2002). The raw event files 
were first examined and corrected for any time anomalies 
(see R02), and then reduced using the XMM-Newton Sci- 
ence Analysis Software v5.3.3. The events from timing mode 
data were extracted in a way similar as in R02. For datasets 
taken in full-frame and small window mode, events were ex- 
tracted from an aperture of ~ 40" centered on the source 
and a background, source-free region; the background sub- 
tracted spectra were obtained by scaling appropriately for 
the size of the two regions. Finally, the event arrival times 
were corrected to the solar system barycenter and to the 
center of the binary system, using the ephemeris of Still et 
al. (2001). 

2.2 Determining the beat phase 

The beat period of Her X-l is known to vary by up to a 
few days (Baykal et al. 1993). In order to determine how 
the epochs of the XMM-Newton observations relate to the 
precession cycle, we extracted the RXTE ASM (2-10) keV 
quick-look light curve from the MIT web site and we deter- 
mined the 35 day phase by computing the time difference 
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Observation 


XMM- 


Exp 


MOS 


Orbital 


35 d 


Date 


Newton 


(ksec) 


Mean 


Phase 


Phase 




Orbit 




Ct/s 






2001 Jan 26 


207 


10 


97.5 


0.20-0.26 


0.17 


2001 Mar 04 


226 


11 


3.6 


0.47-0.56 


0.26 


2001 Mar 17 


232 


11 


11.3 


0.52-0.60 


0.60 


2002 Feb 25 


405 


10.3 


2.5 


0.55-0.62 


0.46 


2002 Feb 27 


406 


11.7 


3.2 


0.71-0.79 


0.51 


2002 Mar 07 


410 


8.8 


5.2 


0.40-0.46 


0.74 


2002 Mar 09 


411 


7.3 


3.4 


0.57-0.62 


0.79 


2002 Mar 15 


414 


7.3 


2.2 


0.09-0.14 


0.96 


2002 Mar 17 


415 


7.3 


130 


0.26-0.31 


0.02 


2002 Mar 27 


420 


4.4 


3.4 


0.17-0.20 


0.31 


2002 Sep 20 


509 


5.9 


2.9 


0.56-0.60 


0.35 


2002 Sep 22 


510 


7.3 


1.8 


0.72-0.77 


0.41 


2002 Oct 04 


516 


10.3 


2.1 


0.78-0.85 


0.77 


2003 Mar 26 


603 


15.6 


2.0 


0.81-0.88 


0.72 


2003 Mar 28 


604 


6.1 


0.5 


0.98-0.03 


0.75 



Table 1. Summary of the XMM-Newton observations of Her X-l. 
For the orbital phasing we use the ephemeris of Still ct al. (2001); 
for the 35 d phasing we define phase 0.0 at the main high state 
turn-on (see text for a discussion). The count rate refers to the 
rate in either the MOS1 or MOS2 detector depending on the 
source intensity. Analysis of the first three datasets has been pre- 
sented in details in R02 and Jimenez-Garatc ct al. (2002). 



between each observation and the start of the previous rise 
to maximum. We then took the period to be the time till the 
next rise to maximum. We estimate the uncertainty on the 
phase of the beat period to be ~0. 01-0. 02 cycles. The result- 
ing beat phase is reported in Table 1, along with the dates of 
the observations and the mean X-ray brightness. Note that 
the last observation has been taken during an eclipse of the 
primary by the secondary star: this corresponds to the point 
with the lowest X-ray count rate in Figure 1. 

The phase coverage of the XMM-Newton observations 
over the 35 day period is shown in Figure 1. For comparison, 
we also show the ASM light curve recorded during a typical 
cycle. We note that, since the actual period of the X-ray 
modulation varies around its nominal value of 35 d, the de- 
termination of the state of the system (either low, main on 
or short on) based on a direct comparison with the ASM 
lightcurve in the top panel must be taken only as indicative. 

As we can see, the beat cycle is reasonably well sampled, 
with most of the exposures taken outside the main-on and 
short-on phases. Three observations have been performed 
during states of high intensity: in addition to the main-on 
and short-on datasets discussed by R02 ($35 = 0.17 and 
0.60) we have a further exposure at $35 = 0.02. We note 
that the ASM lightcurve 




Figure 1. Top Panel: a typical RXTE/ASM (2-10) keV light 
curve of Her X-l over the 35 day beat period. Bottom Panel: 
the mean count rate of Her X-l registered by the EPIC MOS 
detectors. The lowest count rate corresponds to an eclipse of the 
neutron star by the secondary. 



is on average small (less than 10% as estimated by Boroson 
et al. 2001). However, it is modulated over the binary orbital 
period. The UV count rate detected using the two OM filters 
is shown in Figure 2, as a function of the orbital phase. Note 
that observations were not always performed in both filters, 
and in the last two exposures no UV measurements were 
made. The illuminated secondary star is expected to domi- 
nate around ^binary ~ 0.5. As we can see, near this orbital 
phase the UV flux measured by OM reaches a broad max- 
imum, similar to that found using IUE in the 1750-1850A 
band (e.g. Vrtilek & Cheng 1996). However, there is a flux 
decrease at ^binary ~0.5, possibly due to the fact that part 
of the secondary star is shadowed by the accretion disk or 
to a warp in the disk itself. 

The rise from the eclipse {4>binary ~0. 1-0.3) is steeper in 
UVW1 (longer wavelength) than in UVW2. Since, at these 
orbital phases, the accretion disk mainly contributes to the 
UV flux, this is consistent with a scenario in which, after 
the eclipse, the regions of the disk farthest from the neutron 
star come into view first. 



3 THE UV DATA 

The UV emission probably originates from a blend of at 
least two different components: the illuminated face of the 
secondary star and the accretion disk. Their relative contri- 
bution varies during the orbital motion, and in general is 
difficult to disentangle. Unlike the X-rays, the UV emission 
does not show an obvious correlation with the beat period. 
This is in agreement with the fact that the disk contribution 



4 TIMING ANALYSIS 

4.1 Search for modulation close to the spin period 

The spin period of Her X-l has been found to vary from 
~1. 23772 s to 1.23782 s, with alternating phases of spin-up 
and spin-down (Parmar et al. 1999, Oosterbroek at al. 2001). 
Furthermore, it can vary significantly on relatively short 
time scales. The values of the spin period measured using 
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Figure 2. Top Panel: The UVW1 (top panel) and UVW2 (bot- 
tom panel) data folded on the orbital period, data have been 
phased by using the cphcmcris of Still et al. (2001). 



the first three XMM-Newton observations were reported by 
R02: using the PN detector, we found 1.237774 s, 1.237753 s 
and 1.237751 s, respectively. The first measure confirmed 
the slow-down trend monitored by BeppoSax (Oosterbroek 
ct al. 2000) and Chandra (Burwitz, private communication) ; 
although we observed a slight decrease in the two following 
epochs. 

We performed a search for pulsations using a Discrete 
Fourier Transform for all the new datasets, finding a strong 
evidence of pulsations only in one case (revolution 415). 
This is the only new exposure taken outside the low state, 
precisely during the main-on at $ 35 = 0.02. The corre- 
sponding spin period, as obtained from the EPIC PN data, 
is 1.237777(1) s, where the number in parenthesis repre- 
sents the error on the last digit as determined by using 
the Cash (1979) statistic. The spin period is just slightly 
increased with respect to the value measured by XMM- 
Newton in 2001. 

All the remaining datasets showed amplitude spectra 
with only weak or no significant peaks close to the expected 
spin period. In order to investigate further the possible pres- 
ence of a modulation, we folded all the datasets on the most 
recent measured period, i.e. 1.237777 s, by using 25 phase 
bins. We then fitted a constant value to the resulting light 
curve and determined the \ 2 . Therefore, we expect a poorer 
fit in datasets with a stronger spin modulation. Since the 
profile of the folded light curve varies considerably in the 
soft and hard energy band, we split the events into two en- 
ergy ranges: 0.3-0.7 keV and 2-10 keV. The resulting \ 2 is 
shown in Figure 3, as a function of the 35 day phase. For 1 
degree of freedom (the value of the period), the difference 
to the fit is Ax 2 =2.71 at the 90 percent confidence level. 
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Figure 3. Top panel: the mean ASM X-ray light curve over the 
35 day cycle. Middle panel: the x 2 obtained by fitting a constant 
value to the 0.3-0.7 keV folded light curve. Bottom panel: same 
but for 2-10 keV folded light curve. In folding the data we folded 
the data on the spin period determined using a Discrete Fourier 
Transform (main-on or short-on observations) or, when unavail- 
able, the more recent measure of 1.237777 s. The x 2 is computed 
by using 25 phase bins. 



Using this method we find that, outside the main-on and 
short-on states there is evidence for significant modulation 
above 2 keV in the two datasets taken at $35=0.26 and 0.31. 
In the first case (which is one of those already reported by 
R02) the x 2 increases above the 90 percent confidence level 
if we reduce the number of phase bins from 25 to 20. 

Defining the amplitude as (max-min/mean) of a fitted 
sinusoid, we find an amplitude of 11 and 13 percent for 
$35 =0.26 and 0.31, respectively. Previous observations of 
Her X-l during the low state include Coburn et al. (2000), 
who found a spin modulation of ~13 percent (assuming the 
same definition as ours) in the 3-18 keV range using RXTE 
data in an anomalous low state, while Mihara et al. (1991) 
determined an upper limit for the pulsed fraction of 2.4 per- 
cent in the 1.2-37 keV energy range using Ginga data. 



4.2 Spin-resolved light curves 

In Figure 4, we show the spin profiles of Her X-l as derived 
for those observations in which we detected a significant 
modulation in the 0.3-0.7 keV or in the 2-10 keV energy 
band. Due to the uncertainty in the spin period at each 
epoch, the spin phases are not on the same absolute scale. 
Therefore, only the relative phasing between the light curves 
in different energy bands taken during the same observation 
is physically meaningful (see § 4.3). 

The change of the spin profile of Her X-l in different 
energy bands has been well documented in the literature. In 
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Figure 4. Upper panel: the EPIC PN spin profile in the 2-10 keV band as a function of the 35 day cycle. Lower panel: as above but for 
the 0.3-0.7 keV band. In all panels, because of the uncertainty in the spin period in each epoch, the spin phases are not on an absolute 
scale. The two vertical lines at $35 = 0.17 mark the second atypical peak in the interpulsc (hard band) and the small "notch" like feature 
(soft band) discussed in the text. 



order to compare the spin profiles obtained by using XMM- 
Newton with the previously published data, we take as ref- 
erence, in the range above ~ 2 keV, the work of Deeter et 
al. (1998) and Scott et al. (2000). These authors reported the 
variation of the hard X-rays light curve obtained by Ginga 
and RXTE over a wide range of beat phases. 

As we can see, the XMM-Newton light curves obtained 
shortly after the main-on turn on ($35 =0.02) in the 2-10 keV 
energy band closely resemble most of the main-on states pre- 
viously reported. The only exception is that the small peak 
observed during the inter-pulse (at 4> sp in~0.9 in Figure 4) 
appears slightly more sharp than previously reported. At 
the next main-on beat period, $35=0.17, we detect an atyp- 
ical profile, with a second large peak in the interpulse on 
the leading side of the main pulse (the feature occurs at 
(f) sp in~0.05 in Figure 4). The peak becomes more prominent 
as the high state progresses, and has been observed only oc- 
casionally in Her X-l: it is present in the Ginga data taken 
during the main-on state "D" by Deeter et al. (1998), at 
the two beat phases $35 =0.073 and 0.162. However, it was 
not present in the other of their main-on observations nor in 
those discussed by Scott et al. (2000). Deeter et al. (1998) 
noted that such atypical profile is similar to a previous one 
recorded by using Kvant (Sunyaev et al. 1988) and suggest 
it may be due to a reduction in the accretion rate. 

Moving later in beat phase, the detected amplitude of 
the hard X-ray modulation decreases sharply at $35=0.26 
and $35=0.31: here the modulation is broad and quasi- 



sinusoidal. There is evidence of a double peaked structure 
at $35=0.26: R02 show the hard X-ray modulation resolved 
into finer energy bands and find that at this beat phase mod- 
ulation at the hardest energies is broad and single-peaked. 
Note that the profile obtained at $35=0.31 is less certain, 
since it has been obtained by folding data over the most re- 
cent measured period and not using the period measured at 
the time of the observation. However, we have explored the 
effect of folding the data on slightly different spin periods, 
centered around 1.237777 s: the resulting folded light curves 
are very similar. The profile detected during the short-on 
($35=0.60) is similar to previous observations at a similar 
beat phase, although the high signal to noise ratio of the 
XMM-Newton data revealed a more complex structure dur- 
ing the rise of the main peak. 



The high throughput of XMM-Newton allows us to ob- 
tain a light curve with very high signal to noise ratio also in 
the soft X-rays (0.3—1 keV). Consistent with past measure- 
ments, the spin profiles of the soft X-ray band in Figure 4 
show a basically sinusoidal shape. However, in both main-on 
observations the broad peak is cut by a dip which has not be 
seen before. Further, data taken at $ 35 =0.17. show a small 
"notch" like feature preceding the maximum (at sp in=O.4 
in Figure 4). During the short-on ($35=0.60) the modulation 
is broader, with evidence for a fairly complex substructure 
during the leading edge of the peak. 
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Figure 5. The cross correlations for the soft (0.3-0.7 keV) and 
hard (2-10 keV) light curves at four different 35 day phases. 



4.3 The phase shift between the soft and hard 
light curves 

R02 found that the phase shift between the soft and hard 
X-ray light curves, folded on the spin period, varied over 
the beat period. To determine how the dataset taken at 
$35 =0.02 compared with these previous observations we 
cross-correlated the hard and soft light curves (see R02 for 
details). This is the only new observation that can be used 
for such purpose, since the others do not show a significant 
modulation in both soft and hard X-rays. As we can see, 
the two observations taken during the main-on have a very 
similar correlation, with data taken at $35=0.02 showing a 
strong anti-correlation and a phase shift between the max- 
ima of the soft and hard light curves of ~ 130°. This dif- 
ference in phase is slightly smaller than what was observed 
at $35=0.17, and consistent with the general trend over the 
beat period. 



SPECTRAL ANALYSIS: 
OF FE COMPLEX 



THE VARIATION 



5.1 Line variation over the beat cycle 

The integrated X-ray spectra of Her X-l over the 35 day 
cycle have been modeled by many authors (e.g. Mihara et 
al. 1991, Oosterbroek at al. 1997, Coburn et al. 2000, Oost- 
erbroek at al. 2000, Choi et al. 1994, Leahy 2001). In ac- 
cordance with these studies, the main on, short on and low 
states all require a two component continuum with a ba- 
sic emission spectrum and an absorbed component whose 
strength is maximum during the low states. The natural in- 
terpretation is that the 35 day variation is due to the occul- 
tation of the neutron star beam by the intervening accretion 
disk. Model fitting of the XMM-Newton spectra obtained us- 
ing the first three datasets has been reported by R02 and 
confirm this picture. Here, we used the same model for fit- 
ting the newest datasets: since the overall fits do not show 
new characteristics we do not discuss further the details of 
the overall spectral shape. 

Instead, we focus on the behaviour of the X-ray spectra 
near the Fe line component(s), detected in the past around 
~ 6.6 keV. As shown by R02, the central energy of the flu- 
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Figure 6. Variation of the Ka and Fe XXVI line parameters 
along the beat cycle, as measured with EPIC PN. From the top: 
1) mean ASM light curve; 2)-5) central energy, equivalent width 
(EW), width and normalization of the prominent Ka Fe emission 
line; 6)-8) central energy, EW and normalization of the Fe line 
at ~ 7 keV. For those datasets where the second feature is unde- 
tected, we report upper limits to the EW and normalization, at 
the 90 percent confidence level ("v" symbols). 



orescent line varies over the beat period, increasing during 
the main on. However, only three epochs were reported in 
that paper. 

In order to fit the spectra, we restricted the energy range 
to 4-10 keV and fitted an absorbed power law with one or 
more Gaussians. The power law and Gaussian parameters 
were then fixed and the energy range 0.3-4 keV included. 
The absorption component was allowed to vary and we in- 
cluded a blackbody component and a neutral partial cover- 
ing model if required. This allows us to fix the absorption 
components at an appropriate level, which is crucial since, 
for large values of the column density, a characteristic ab- 
sorption feature appears at ~ 7.1 keV. This feature may 
affect the identification of the emission Fe lines, therefore it 
has to be modeled carefully. 

Consistent with R02, we found that a prominent fluo- 
rescent line is detected at all beat phases (apart from obser- 
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vation taken during at mid- eclipse ) , with a central energy 
that varies between ^6.4-6.6 keV and is highest during the 
main-on. However, the most striking result is the detection 
of a second component, an Fe XXVI line at ~ 7 keV, which 
is only present outside the main-on. This line was not re- 
ported by R02, but a closer analysis of the two observations 
taken at $35 = 0.26 and 0.60 revealed its detection in both 
datasets. The line parameters are reported, as a function of 
the beat phase, in Figure 6. 

We find that the equivalent width (EW) of the Ka Fe 
line shows a variation through the beat cycle, with no trend 
being visible. This is likely to be due to a varying contin- 
uum, since the Ka Fe line normalisation is highest during 
the main-on state while it remains low and almost constant 
at other 35 day phases. Similarly, the line width is greatest 
during the main-on state. The energy of the second com- 
ponent remains approximately constant at ~7.0 keV, ex- 
cept during the short-on state where it is significantly lower 
(~6.85 keV). The feature is strongest (larger EW) during 
the states of low X-ray intensity, although it is undetectable 
in some of the exposures in this phase range. On the other 
hand, its normalisation is strongest during the short-on. In 
some of the observations, we do not find evidence for a sig- 
nificant Fe XXVI line. In order to determine upper limits 
to EW and line normalisation, we added a Gaussian com- 
ponent to the best fit spectral model, by fixing the central 
energy at 7 keV and the width at 0.1 keV. The upper limits 
obtained for the datasets taken during the main-on are con- 
sistent with the line being undetected because of the strong 
continuum emission. 

In order to show the significance of both the Fe features, 
we show in Figure 7 the X-ray spectrum between 5-9 keV 
for $35=0.02 and $35 =0.79. The solid line is the best fit 
spectral model with the normalisations of the two Gaussian 
components set to zero. We find no evidence for the presence 
of a Compton downscattcred shoulder (expected near Eq — 
Eo/(m e c + -Bo) where Eq is the centroid of the line) of the 
kind recently observed in HETGS Chandra spectra of the 
high mass X-ray binary GX 301-2 (Watanabe et al. 2003). 



5.2 Pulse-resolved spectroscopy 

The variation of the Fe Ka line over the neutron star spin 
period and its correlation with the thermal or non-thermal 
spectral components was studied by R02. We found no ev- 
idence for a significant variation in the line parameters at 
$35 =0.26 or 0.60. In contrast, the observation made dur- 
ing the main-on ($35=0.17) showed a clear correlation be- 
tween the EW of the 6.4keV line and the soft X-ray light 
curve. This supported the idea of a common origin for the 
~ 6.4 keV fluorescent line and the thermal flux. 

The new data taken during the main-on ($35 =0.02) can 
be used to test this scenario further. We phased all the events 
on the spin period determined in §4 and extracted spin-phase 
resolved spectra. Details of the model fitting are as in §5.1. 
Results are reported in Figure 8, together with the soft and 
hard X-ray light curves. 

The shape of the EW variation is broadly similar to the 
profile at $35=0.17 (R02), although that observation shows 
a more 'top-hat' profile. Consistent with R02, we find in the 
main-on state that both the line normalization and the EW 
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Figure 7. The spectral region around 6.6keV from <3?35=0.02 
(top) and <I>35 =0.79 (bottom). The y-axis is the normalised 
counts, i.e. the number of counts per energy bin per second, as 
measured by EPIC PN. In both panel a solid line shows the best 
fit after the normalisation of the one or two Gaussian components 
have been set to zero. 



follow the soft X-ray light curve and are anti-correlated with 
the hard, non-thermal flux emitted above 2 keV. 



5.3 Variation of Fe complex over the orbital 
period 

Further information about the origin of the Fe Ka 6.4keV 
line can be obtained by observing its variation over the or- 
bital period. This is shown in Figure 9, for all observations 
taken outside the main-on. The average line normalization 
is 2 — 4 x 10~ 4 counts/cm 2 /s, with larger peaks at ~ 10~ 3 
counts/cm 2 /s and a dip at 4>binar y = 0.43 — 0.56. (The line 
normalisation during the main-on state is greater by a factor 
if >5-10 compared to the low state). The line is not detected 
during the middle of the eclipse, at which point the upper 
limit on the line normalisation is 3.7 x 10 -5 counts/cm 2 /s. 

In the previous section we showed that during the main- 
on state the line normalisation phased on the spin period 
is correlated with the soft X-ray light curve suggesting a 
common origin for the thermal soft X-ray emission and the 
6.4keV line. We now find that in the low state the line nor- 
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Figure 8. The observation made using the EPIC PN detector at 
$35 =0.02. From the top panel: the light curve in the 0.3-0.7 kcV 
band; the 2-6 keV band; the Fe Ka line normalisation and EW 
as a function of the spin phase. 



malisation is correlated with the orbital period. We go onto 
explore this further in §6.3. 



6 DISCUSSION 

We have reported the results based on EPIC and OM data 
for a set of XMM-Newton observations of Her X-l. Con- 
sistent with past reports, we find that the spin modulation 
of the neutron star is most prominent during the main-on 
and short-on state, becoming not significant during the low 
state. These data sets allow us to follow the change in the 
relative phasing of thermal and non-thermal emission and 
the change in the pulse profile over the beat period. The 
large collective area of XMM-Newton and its good response 
down to soft X-rays also revealed a complex substructure in 
the light curve in the softest and less well studied energy 
band. 

Data taken during the main on, at $35 = 0.02, show a 
common phase modulation between the fluorescent Fe Ka 
line and the soft emission, confirming the results found by 
R02 at $35 = 0.17. Moreover, for the very first time in CCD 
non grating spectra, we report evidence for a second Fe line, 
with central energy near ~ 7 keV, which is only visible out- 
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Figure 9. The Fe 6.4keV Koj line normalisation as a function 
of the orbital phase, for all observations taken outside the main- 
on. Data have been taken with EPIC PN. Only an upper limit 
has been obtained from the observation taken during the eclipse 
(arrow) . 



side the main-on state. We now proceed by discussing some 
physical interpretations for our results. 



6.1 The spin-resolved light curves: evidence for 
structure in soft X-rays (0.3 — 1 keV) 

The complex, multipeaked structure of the pulse profile of 
Her X-l above ~ 1 keV and its evolution pattern over the 
beat cycle is the subject of several studies. The basic fea- 
tures of the pulse profile are nominally explained by the 
oblique rotator model (Lamb, Pethick & Pines 1973), but 
there are still many unresolved questions regarding the de- 
tails of the mechanism, such as whether the radiation pat- 
tern is best described by a superposition of 'fan' and 'pen- 
cil' beams and whether the obscuring material that modi- 
fies the visibility of the direct beams from the neutron star 
contributes to the pulse formation. Particularly successful 
in explaining the pattern observed by RXTE and Ginga 
is the accretion column model by Scott et al. (2000). This 
model invokes the successive obscuration of a direct beam 
(that originates from the polar caps of the neutron star), 
two fan beams focussed in the antipodal direction (possibly 
due to cyclotron backscattering) and a more extended scat- 
tered halo. As the main-on progresses, the inner edge of the 
accretion disk cover first one of the fan beam components, 
then the direct pencil beam and then the second fan beam 
component. 

In contrast, the pulsation detected below ~ 1 keV us- 
ing EXOSAT (Ogelmen & Triimper 1988), ROSAT (Mavro- 
matakis 1993), and BeppoSax (Oosterbroek et al. 1997, 
2000, 2001) is broad and quasi-sinusoidal. This emission is 
thermal and probably originates in the inner edge of the 
disk, that partially intercepts and reprocesses the hard X- 
rays from the neutron star. XMM-Newton data taken during 
the main-on have an excellent signal to noise ratio and allow 
us to resolve a substructure in the soft X-ray light curve. The 
0.3-0.7 keV pulse profile observed at $35=0.02 clearly shows 
the presence of two separate pulses, at 4> sp in ~ 0.0 — 0.1 and 
4>a P in ~ 0.3 — 0.4 in Figure 4 (the value of the absolute 
spin phase is arbitrary) . Similar features are also evident in 
data taken later in the beat cycle: at $35=0.17 the pulse 
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profile shows two maxima (at sp i„ ~ 0.7 and <p 3 pi„ ~ 0.9 
in Figure 4) and a small "notch" like feature preceding the 
main pulse (at 4>apin=0A in Figure 4). During the short on 
($35 =0.60), the soft X-ray modulation is broader, more si- 
nusoidal, and exhibits a complex substructure in the leading 
edge of the pulse followed by a faster decay. 

Whether these features reflect an intrinsic complexity in 
the thermal and geometrical properties of the reprocessing 
region or keep memory of the structure of the illuminating 
beam(s) is unclear. However, it is worth noticing that during 
all three beat phases with bright emission, the features de- 
tected in the soft light curves qualitatively resemble those of 
the corresponding hard band. The two maxima observed in 
the softer band during the main-on ($35=0.02, 0.17) have 
a phase separation similar to that of the two main peaks 
observed in the 2-10 keV range, i.e. A(f> spin ~ 0.2 — 0.3. 
Similarly, the small "notch" preceeding the main pulse at 
$35=0.17 may be associated with the atypical pulse ob- 
served at higher energies (both precede the main peak by 
Atfispin ~ 0.25 — 0.3). If these associations are correct, i.e. if 
at least some of the main characteristics of the neutron star 
beam are not "washed out' during the reprocessing of hard 
X-rays, this means that the radiative time scale in the inner 
edge of the disk is much smaller than the photon travel time. 
Moreover, the reprocessing region cannot be too extended in 
size. This may be explained by a high inclination angle of 
the disk and/or a highly warped inner edge, in which case 
only a fraction of the inner region of the disk intercepts the 
neutron star beam. Alternatively, the reprocessing of the 
neutron star beam may be dominated by a shocked, opti- 
cally thick hot spot instead of being distributed over the 
entire inner edge. A similar situation is observed in SW Sex 
stars (see e.g. Dhillon et al. 1997, Groot et al. 2000 and 
their Figure 6) , where the hot spot region is formed through 
a shock that occurs along the stream trajectory, but well 
inside the accretion disk. 

6.2 The pulse period and the change in the phase 
shift between soft and hard emission 

Given the complexity of the source, pulse-phase spec- 
troscopy is of paramount importance in disentangling the 
different spectral components observed in Her X-l. Past ob- 
servations using Einstein (McGray et al. 1982) and Bep- 
poSax (Oosterbroek at al. 1997, 2000) showed that, during 
the main-on state, the maximum of the thermal component 
and the power- law components are shifted by ~ 250°. The 
departure from a phase shift of 180° is usually associated 
with a disk having a tilt angle, under the assumption that 
the soft X-rays originate in the layers of the disk that inter- 
cept and reprocess the neutron star beam. The first evidence 
for a change in the pulse difference, as measured at three dif- 
ferent $35 using XMM-Newton, has been reported by R02. 
R02 found that during the main-on state ($35=0.17) the 
soft and the hard X-ray pulse profiles were strongly anti- 
correlated and that the phase shift between the soft and 
hard maxima was ~ 150°. This was in contrast to previous 
observations made during the main-on. Moreover, the phase 
shift between these components was significantly different in 
the XMM-Newton observations performed at $35=0.26 and 
0.60. This led to the suggestion that, during the first three 
exposures, the tilt angle of the accretion disk was chang- 



ing substantially (R02). In contrast with past observations 
(Oosterbroek et al. 2000), R02 did not detect a symmetric 
state at $35 ~ 0, 0.5. If the reprocessing region is located at 
or near the inner region of the disk, such lack of symmetry 
may indicate that the disk is strongly warped (Heemskerk 
& van Paradijs 1989). In order to further monitor the phase 
shift over the beat cycle, we presented a new measure of the 
pulse difference, taken during the main-on ($35 = 0.02). We 
find that the trend detected by R02 is confirmed, and the 
change in the tilt angle appears to be smooth and continuous 
over the beat cycle. 

Possible recent changes in the structure of the inner 
layers of the disk have been also discussed by Oosterbroek 
et al. (2001), who detected a rapid spin down in Her X-l 
following an anomalous low state. This spin down was not 
accompanied by a significant simultaneous change in the 
mass accretion rate (Oosterbroek et al. 2001). Instead, it 
was more plausibly associated with changes in the structure 
of the outer magnetosphere of the neutron star and/or in the 
inner disk layers. These observations impose important con- 
straints on any theory modeling the accretion torques and 
phenomena related to the inner edge of the disk in Her X-l. 

6.3 The Fe Ka line 

The strong emission line at ~ 6.4 keV is detected in all 
XMM-Newton observations, with larger broadening and nor- 
malization during the main-on. Past reports of Fe features 
in the spectra of Her X-l include an ASCA exposure (Endo 
et al. 2000) and a set of 63 Ginga observations analyzed 
by Leahy et al. (2001). The good coverage of the beat pe- 
riod allowed Leahy et al. (2001) to undertake a systematic 
study of the line variation along the precession cycle. They 
found a broad feature near ~ 6.5 keV, with the EW vary- 
ing considerably between main on, short on and low state. 
The mean values for each state are 0.48 keV, 0.66 keV, and 
~ 1.37 keV, respectively. Similar correlations are not evident 
in the XMM-Newton data, where instead an additional Fe 
XXVI line, possibly a coronal signature, is detected during 
the low states (see § 6.4). 

Nevertheless, we find evidence for a significant variation 
of the line energy over the 35 day period: the Fe line emission 
originates in near neutral Fe (Fe XIV or colder) in the low 
and short-on state observations, whereas in the main-on the 
observed Fe Ka centroid energies (6.65 ±0.1 keV and 6.50 ± 
0.02 keV as measured using PN at $35 = 0.02 and 0.17) 
correspond to Fe XX- Fe XXI (Palmeri et al. 2003). The line 
centroids observed using the EPIC PN deviate by 4<r from 
the 6.40 keV neutral value. This has been already noticed 
by R02, who suggested three possible explanations for both 
the line broadening and the centroid displacement: 

(i) an array of Fe Ka fluorescence lines exists for a variety 
of charge states of Fe (anything from Fe I-Fe XIII to Fe 
XXIII); 

(ii) Comptonization from a hot corona with a significant 
optical depth for a narrower range of charge states centered 
around Fe XX. Similar line broadening has been observed in 
some low mass X-ray binaries with ASCA (Asai et al. 2000), 
and Her X-l spectra exhibit one of the larger equivalent 
widths within the class; 

(iii) Keplerian motion: if this is responsable for the line 
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broadening, the Keplerian velocity measured at $35 = 0.02 
and 0.17 is ~ 15500 and ~ 13000 km/sec, respectively. This 
gives a radial distance of ~ 2 — 3 x 10 8 cm (for a neutron star 
of 1.4 solar masses), which is close to the magnctospheric 
radius for a magnetic field of ~ 10 12 G. 

In addition to these possibilities, we also notice that 
the region responsible for the Fe Ka line emission is likely 
to be different for lines observed at different beat phases. 
The Fe lines have different energies, flux, and broadening in 
the main-on compared to the low state. While data taken 
during the main on clearly suggest a correlation between the 
fluorescent Fe Ka line and the soft X-ray emission, the same 
is not explicitly evident in data taken during the low state. 
At such phases, instead, the line is a factor > 5 weaker and 
is clearly modulated with the orbital period. 

The correlation between the UV and the neutral Fe K 
point to a common origin with the fast rise UVW1 emission. 
Still, it is hard to distinguish the contributions from the 
companion and the disk to the Fe K line emission. A large 
contribution from the companion would explain why the Fe 
K line is so strong in the low state of Her X-l compared to 
other accretion disk sources: in Her X-l the companion is 
much larger (2.3 Mq) than for other LMXBs. However, due 
to the fast rise of the Fe K flux with the orbital phase, the 
disk origin scenario is still strong. 

A possible contribution to the Fe Ka emission may 
arise from relatively cold material in an accretion disk wind, 
such as commonly observed in cataclysmic variables (see e.g. 
Drew 1997). Although no consensus has been reached in the 
case of low-mass X-ray binaries, in the case of Her X-l evi- 
dence of outflowing gas, possibly a wind or material ablated 
from HZ Her, has been reported by Anderson et al. (1996) 
and Boroson et al. (2001). UV lines observed with the FOS 
and STIS spectrographs on the Hubble Space Telescope per- 
sist even in the middle of the eclipse, when the X-ray heated 
atmosphere of the normal star and the accretion disk should 
be entirely hidden from view. The velocities inferred from 
the broadening of the N V lines are not constant over a time 
scale of ~ 1 yr, while a comparison between observations 
taken at (pbinary =0.10, 0.21 and 0.29 show they are stable 
over the orbital phase. Between these phases, the velocity of 
the neutron star is expected to change by « 60 km/s. In this 
respect we note that the Fe Ka line was not been detected 
by XMM-Newton during the middle of the eclipse, and the 
upper limit on the line flux is ~ 10 or less of that measured 
outside the eclipse. This in turn sets an upper limit on a 
potential contribution to the emission of this line by a wind 
or some sort of circumstellar material. There is no Doppler 
signature of a wind in the HETG spectrum of the Fe K line 
(Jimenez-Garate et al. 2003), and the wind would have to 
be enclosed by the Roche lobe due to the absence of Fe K 
in eclipse. 

On the other hand, taken as a body the data reported 
here suggest a complex origin for the overall emission of the 
Fe Ka line. To our knowledge, a complex of lines which in- 
clude all ionization states from Fe XVIII to XXIV Fe Ka has 
not been observed in any astrophysical source. This may still 
indicate an outflow of relatively cold gas or some complex 
dynamics in the disk/magnetosphere interface. Such phe- 
nomena should be time-dependent and may be monitored 
in the future using Astro-E2. 



6.4 The ~ 7 keV Fe line 

Our observations have revealed an Fe XXVI line at ~ 7 keV, 
detected during low-state. The line is not detected during 
the main-on when the X-ray spectrum is dominated by the 
strong continuum emission. This feature has not been de- 
tected by Leahy et al. (2001) using Ginga data. However, 
Ginga has a relatively low spectral resolution around these 
energies (18 percent at 6 keV). ASCA observed the source 
with better spectral resolution (2 percent at 5.9 keV using 
the SIS detectors) during an early main on state, and based 
on these data Endo et al. (2000) reported the first resolu- 
tion of the broad feature into two components, at 6.41 keV 
and 6.73 keV. However, this double structure has not been 
detected in other observations taken during the main-on nor 
is it present in XMM-Newton data at similar phases. 

The EPIC detectors on board XMM-Newton have a 
spectral resolution similar to ASCA at 6 keV. However, 
XMM-Newton has a much greater effective area allowing 
lines to be detected with a greater confidence, in particular 
during the states of low intensity of the binary system. The 
feature at ~ 7 keV has been observed by XMM-Newton 
for the first time over several 35 day phases. Also, it has 
been confirmed by a Chandra HETGS observation of the 
source (the only one made during the low state) taken at 
$ 35 = 0.44 - 0.46 (Jimenez-Garate et al. 2003). The fea- 
ture cannot be produced by fluorescence, and it is more 
likely to be a Fe XXVI line originating in widely extended 
photo-ionized plasma. On the other hand, RGS data taken 
during the low and short-on states also show the presence 
of photo-ionized gas (Jimenez-Garate et al. 2002). Grating 
spectra exhibit several narrow recombination emission lines, 
the most prominent being C VI, N VI, N VII, O VII, O VIII 
and Ne IX. The line ratio G = (/ + i)/r, as computed for 
all the helium-like ion complexes, is G ~ 4, which indicates 
that photoionization is the dominant mechanism.* More- 
over, RGS spectra shows two radiative recombination con- 
tinua of O VII and N VII, consistent with a low temperature 
of the emitting plasma (30000 K < T < 60000 K, Jimenez- 
Garate et al. 2002), which furtherly validate the photoion- 
ization model. None of these features is detected during the 
main-on state (where instead there might be evidence for a 
broad O VII i line, and perhaps O VIII Ly alpha). The re- 
combination X-ray line emission are not likely to originate in 
HZ Her, due to the absence of UV induced photoexcitation 
signatures in the He-like triplets (observed with HETGS, 
Jimenez-Garate private comm.). 

Instead, as suggested by Jimenez-Garate et al. (2002), 
an extended, photoionized accretion disk atmosphere and 
corona may be responsable for such features. Jimenez- 
Garate et al. (2001) computed the spectra of the atmo- 
spheric layers of a Shakura-Sunyaev accretion disk, illumi- 
nated by a central X-ray continuum. They found that, un- 
der these conditions, the disk develops both an extended 
corona which is kept hot at the Compton temperature, and 
a more compact, colder, X-ray recombination-emitting at- 
mospheric layer (see Figure 9 in Jimenez-Garate et al. 2002). 
The parameters of the compact atmospheric layer are consis- 

* Here i, r and / denote the intercombination line blend, the res- 
onance line and the forbidden line; for collisionally ionized gases 
it is G < 1 (Licdahl 1999, Porquet & Dubau 2000). 
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Figure 10. Electron density plotted versus the size of the region 
enclosing the Fe XXVI line emitting plasma. The three lines are 
the constraints discussed in the text and the shadowed region is 
allowed. 

tent with the constraints to the narrow line emitting region, 
as derived from spectroscopic analysis and modeling of the 
RGS features. 

Interestingly, the Fe XXVI line detected by XMM- 
Newton may be a signature of the hottest external layers of 
the disk corona, which are located above the recombination- 
emitting layers. In order to detect a line at ~ 7 keV, an 
ionization parameter log£ > 3.3 is required (Kallman and 
McCray 1982, their Fig. lb). Using the flux of power law 
component reported by R02 (i.e. 3.65 x 10 _s erg/cm 2 /s at 
$35 = 0.17) and a source distance of 6.6 kpc, the luminosity 
of the direct beam is L x = 1.9 x 10 3S erg/s. This gives 

nr 2 = Ij- < 9.7 x 10 34 cm -1 , (1) 

where n is the particle density and r is the size of the re- 
gion enclosing the Fe XXVI line emitting plasma. A second 
constraint can be derived from the emission measure. Run- 
ning a detailed multi-temperature plasma model in XSTAR 
gives an emission measure of n 2 V > 6.1 x 10 57 cm" 3 for the 
7 keV line, for a flux of 2 x 10~ 4 photons/cm 2 /s (the emis- 
sion measure scales linearly with the line flux, and the latter 
varies between 1 — 2 x 1CF 4 photons/cm 2 , see Figure 6)t 
. The average temperature, computed by using a main-on 
custom ionization spectrum in XSTAR as well as HULLAC 
recombination emissivities (D. Liedahl, private comm.), is 
kT ~ 580 eV. Moreover, if the line emitting plasma is lo- 
cated outside the Alfven shell, r > 3.6 x 10 7 cm. These 
three constraints define an allowed region in the plane n — r, 
which is shown in Figure 10. The density of the hot coronal 
layers in the disk models by Jimenez-Garate et al. (2001) is 
10 15 - 10 16 cm" 3 if the radius is r ~ 10 s - 10 10 cm, therefore 
within the limits shown in the figure. 

In summary, the variability of the Her X-l spectrum 
lends support to the precession of the accretion disk. The 
evidence for the disk identification relies on the modeled 

t The value of the emission measure is consistent with calcu- 
lations by Kallman ct al. (2003), who found 10 57 cm" 3 for an 
accretion disk corona at 7 kpc and a Fe XXVI line with a roughly 
similar flux level. 



structure and spectra from a photoionized disk (Jimenez- 
Garate ct al. 2001), which is in agreement with the limit set 
spectroscopically on the density of the low-energy lines emit- 
ting region. Interestingly, the same model naturally provides 
a candidate for the region emitting the Fe XXVI line and 
again the computed values of the density agree with the con- 
straints inferred from the line parameters. Taken together, 
these facts strengthen the interpretation of the low state 
omission in term of an extended source and open the ex- 
citing possibility to monitor spectroscopically the different 
atmospheric components of the disk during the transition 
from the low to the high state. 
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